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ABSTRACT 



A Naval Postgraduate School spacecraft design class 
proposed a multiple beam, phased array, microstrip antenna as 
part of the preliminary design of a low earth orbit communica- 
tion satellite. The antenna must provide coverage over the 
satellite ' s entire field of view while both uplink and down- 
link operate simultaneously on the same L-band frequency. 

This thesis assesses the feasibility of the antenna 
proposed in that preliminary design. Design tradeoffs for a 
microstrip array constrained by both available surface area 
and a limited mass budget are examined. Two different 
substrate materials are considered in terms of weight and 
I performance . Microstrip patch theory is applied to array 
element design and layout and antenna array theory is applied 
to determine phase and amplitude coefficients. The focus of 
the design is on obtaining the desired beam shape and orienta- 
tion, given antenna size constraints. A corporate feed method 
is discussed and a general design presented. 

Antenna performance is predicted through the use of a 
computer model based on Modal Expansion theory and results are 
plotted in a series of graphs which demonstrate the limita- 
tions of the proposed design. 
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I. INTRODUCTION 



A. BACKGROUND 

In June 1991, Loral Cellular Systems Corporation of Palo 
Alto, California submitted a license request to the Federal 
Communications Commission for authority to build a world-wide, 
mobile communications system [Ref. 1] . Their proposed system, 
called GLOBALSTAR, is a satellite system designed to provide 
Radio-Determination Satellite Services (RDSS) for real time 
position location and tracking, as well as voice and data 
services . The system is intended to be integrated into the 
existing Public Switched Telephone Network (PSTN) and private 
cellular networks in order to provide telephone communications 
to users anywhere in the world. 

In a project sponsored by the NASA/Universities Space 
Research Association Advanced Design Program, students at the 
Naval Postgraduate School, Monterey, CA performed a prelimin- 
ary design [Ref. 2] for a satellite system based on the 
proposed GLOBALSTAR system. The focus of this thesis is on 
the design and performance analysis of a phased array L-band 
antenna, coirposed of microstrip patches, which was proposed as 
part of this preliminary design. The scope is limited to the 
design of the antenna panels and microstrip patches and to the 
performance analysis of overall antenna. While the beam 
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forming network will be mentioned, a detailed analysis of the 
beam forming network is outside the scope of this thesis. 
Similarly, certain mechanical details such as panel hinging 
and antenna deployment mechanisms as well as specific 
fabrication methods will not be considered in detail. 

B. SYSTEM (CONCEPT 

Using this system, a mobile telephone user will have 
access to world wide coverage by accessing the PSTN via 
satellite. Using a special mobile telephone unit, a user 
communicates directly with a satellite via an L-band radio 
link. The satellite receives signals from mobile users and 
retransmits these signals to an earth based gateway via an S- 
band trunk link. In the reverse process, signals from the 
gateways are received by the satellite via the S-band link and 
then retransmitted to mobile users via the L-band link. 
Gateways are distributed across the globe and are connected to 
the PSTN. 

Acquisition, synchronization, and satellite beam hand off 
coordination is provided for the gateways by a Satellite 
Operation Control Center (SOCC) as depicted in Figure 1.1. 
The information provided by the SOCC allows each gateway 
station to track and communicate with the satellites in its 
field of view. Ephemeris data is provided to the SOCC by 
several Tracking Telemetry and Control Stations (TT&C) which 
obtain the data through satellite observations. Overall 
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communications network management is provided by a Network 
Control Center (NCC) . 




Figure 1.1 System Architecture 



C. SYSTEM DESIGN 

An orbital altitude of 750 nm gives a satellite field of 
view of approximately 108 degrees (assuming a minimum 
elevation angle of 10 degrees for the satellite observer) . 
This allows a constellation of 48 identical satellites in 
circular orbits to provide continuous coverage between 75 
degrees North latitude and 75 degrees South latitude. 
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The constellation is divided into eight orbital planes at 
an interval of 30 degrees with six satellites in each plane. 
There is no provision for on-orbit spares. Communication 
between the satellite and gateways is accomplished via an 
uplink and downlink at center frequencies of 6533.25 MHz and 
5207.75 MHz respectively, each with a bandwidth of 16.5 MHz. 
A mobile user communicates with a satellite via an L-band link 
operating at a center frequency of 1618.25 MHz and a bandwidth 
of 16.5 MHz. A summary of the link budget as specified by 
[Ref. 1] is contained in Appendix A. 

Each satellite provides communications coverage over its 
entire field of view through the use of six non -overlapping 
elliptically shaped spot beams. To maximize the amount of 
time a mobile user is illuminated by any one particular beam, 
beams are oriented such that their major axes are parallel to 
the spacecraft velocity vector. Figure 1.2 shows a typical 
coverage pattern for one satellite. 

In order to minimize interference between channels, the 
six beams operate in pairs such that maximum beam separation 
is achieved on the ground. This is accoirplished by switching 
beams on and off in pairs so that opposite pairs are active 
for 10 ms out of a 60 ms duty cycle. During the duty cycle, 
beam pairs alternate between transmit and receive functions. 
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Figure 1.2. L-Band Beam Orientation 
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D. ANTEINNA PERFOPMANCE REQUIREMENTS 

The L-band antenna for this system must not only provide 
adequate gain, but must also correctly shape each of the six 
spot beams while maintaining low sidelobes . The antenna must 
also allow the communications payload to maintain a constant 
bit error rate (BER) . Bit error rate is a function of carrier 
to noise ratio (C/N) and C/N is a function of path length. 
Because antenna beams must be statically scanned off center to 
cover the FOV, the path length for each of the six antenna 
beams is different. Path lengths vary from a minimum of 1748 
km (measured along antenna boresight) for the inner most 
beams, to a maximum of 9066 km for the outer most beams. In 
order to maintain a constant BER, an isoflux antenna is 
required. Table I-l summarizes the L-band antenna performance 
parameters . 



TABLE I-l. ANTENNA PARAMETERS 



fo=1618 .25 
BW=16 . 5 

(MHz) 


Beam 
1 & 6 


Beam 
2 & 5 


Beam 
3 & 4 


Center 


45° 


30° 


H 

O 

o 


BW (-4 dB) 


25° X 110° 


35° X 120° 


40° X 120° 


BW (total) 


40° X 120° 


50° X 120° 


70° X 140° 


Gain (dBi) 


6.5 


5.5 


4.0 



Note: Satellite FOV is 108° 
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II. DESIGN CONSIDERATIONS 



A. SPACECniAFr CONFIGURATION 

1 . Overview 

The spacecraft is designed to be placed into orbit 
using the Delta 7925 launch vehicle. A Satellite Launch 
Dispenser (SLD) was designed so that a single launch vehicle 
can place an entire plane of satellites into orbit at one 
time. The SLD is essentially a chest of drawers where the 
spacecraft are the drawers. The spacecraft slide into the SLD 
along rails and compress a spring when properly seated. Each 
spacecraft is then held in place by explosive bolts . At the 
appropriate time during orbit insertion, the explosive bolts 
are fired and a spacecraft is ejected into the proper orbit. 

2 . Size Constraints 

The Delta payload shroud and the SLD impose serious 
size constraints on overall spacecraft size. In order to fit 
into the space provided by the SLD, the spacecraft dimensions 
can not exceed 1.4 x 1.4 x 0.6 meters. The space allowed for 
the L-band antenna alone is 1.4 x 1.4 x 0.1 meters in the 
stowed configuration. The preliminary design allowed an 
antenna mass budget of 10 kg however this proves to be 
unrealistic . 
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B. ANTENNA CONFIGURATION 



In order to achieve beam steering and beam shaping 
characteristics necessary to meet the antenna performance 
requirements, a planar phased array antenna similar to the L- 
band antenna of [Ref. 3] was considered. A planar array 
allows both beam shaping and beam steering in any direction. 
For simplicity, rectangular microstrip patch array elements 
were chosen. A microstrip antenna is both lightweight and 
easily fabricated using common printed circuit board etching 
techniques [Ref. 4]. Disadvantages usually associated with 
microstrip antennas, specifically large surface area require- 
ments and narrow bandwidth [Ref. 5] , were found to be insigni- 
ficant in this application. Folding array panels are employed 
so that the antenna can be stowed in the space available and 
conveniently deployed via a spring-loaded mechanism once the 
spacecraft is in the proper orbit . Each array panel 
corresponds to one of the six antennas necessary to form the 
appropriate beam . Figure 2 . 1 shows the phased array antenna 
in the deployed configuration. 

C. MICROS'ERIP PATCH THEORY 

The simplest and most common type of microstrip antenna is 
the rectangular patch mounted on a dielectric substrate of 
thickness h, shown in Figure 2.2. 
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1 . Radiation Mechanisms 

The radiation mechanism of a microstrip patch can be 
modeled in one of several ways, each with its advantages and 
disadvantages . 

a. llie Wire Grid Model 

A microstrip patch can be modeled as a grid of 
fine wires [Ref. 6] . The currents on the wires are solved for 
numerically and when found, provide accurate information about 
the antenna pattern and input impedance. The disadvantage of 
this method is that it requires an extraordinary amount of 
computer resource [Ref. 7] . 
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b. The Modal Expansion Model 



The modal expansion model, which is valid for 
electrically thin (h«?i) substrates, represents the antenna as 
a thin TM mode cavity bounded by magnetic walls [Ref. 5] 
within which the dominant TM mode is excited along with other 
non-resonant modes. As long as the dominant mode is suffi- 
ciently excited, this model provides accurate performance 
information [Ref. 7], including the effects of feed location. 
The ability to account for a feed probe and handle complex 
patch shapes are the chief advantages of this model . 
c. The Transmission Line Model 

The transmission line model considers the patch to 
be a line resonator with radiation due to fringing fields at 
the open circuited ends [Ref. 6] . The radiator consists of 
two radiating slots separated by a distance L. The slots are 
formed at the edges of the patch element, between the patch 
and the ground plane below. Field variations along the 
radiating edges are ignored. The chief disadvantage of the 
transmission line model is that it is only valid for rectan- 
gular patches. The input impedance of the transmission line 
model is also highly dependent on feed location and fails to 
accurately predict input impedance for all feed locations 
[Ref. 7] . However, formulas have been developed to predict 
input impedance for any feed location [Ref. 6] . 

Despite its drawbacks, the transmission line 
model's main advantage is, that for rectangular patches, it 
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provides simple design formulas which give reasonably accurate 
performance predictions . 



2 . Patch Parameters 

The dimensions of a rectangular patch, as well as 
bandwidth and gain, are determined by the operating frequency 
of the antenna, the relative dielectric constant, and 
thickness of the substrate material . The following formulas 
are' based on the transmission line model, 
a. Width and Lengrth 

The width and length of a rectangular microstrip 
patch are given by. 



where 

c = speed of light (m/s) , 

= operating frequency (MHz) , 

Cj. = relative dielectric constant, 
Cg = effective dielectric constant. 




(cm) 



( 2 . 1 ) 



L = — - 2Ai (cm) 




( 2 . 2 ) 




( 2 . 3 ) 
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(2.4) 



Ai = 0.412/2 



(e_+0.3) (-^ + 0.264) 
h 

(e^-0.258) (-^+0.8) 
® h 



and 



h = dielectric thickness (cm) [Ref. 6] 



b. Bandwidth 

The bandwidth of a microstrip antenna is defined 
in terms of the antenna ' s quality factor (Q) as follows : 



BW = 



VSWR-1 



(2.5) 



Q/vsm 

where VSWR is less than a specified value (2:1 or 1.5:1 for 
example) and VSWR=l at the operating frequency [Refs. 5 and 
6] , the bandwidth is dependent on both the relative dielectric 
constant, e^, and thickness of the substrate. Thicker sub- 
strates and lower values of give larger bandwidths. 

Typically, microstrip antennas have bandwidths on the order of 
a few percent of the operating frequency, 
c. Gain 

Gain is also affected by substrate thickness and 
relative dielectric constant . Gain is inversely proportional 
to and directly proportional to substrate thickness. To 
achieve the necessairy gain, microstrip patches can be arranged 
to fom linear and planar arrays. 



I 
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D. ARRAYS 



The radiation pattern of a single array element is 
relatively wide and typically has low gain. In order to 
achieve higher directivity and additional gain, antenna 
elements can be arranged to form linear or planar arrays. 

1. Theory 

The total field of an array is given by the vector 
addition of the fields radiated by individual elements [Ref. 
8] . The required directivity is achieved by arranging 
individual element excitations such that fields produced by 
each element interfere constructively in the desired direction 
and interfere destructively in all other directions. The 
radiation pattern of an array is obtained by multiplying the 
field produced by a single element by an array factor (AF) . 
The array factor for an array of N identical elements is given 
by [Ref. 8] as 



N 



AF = J„e 






n=l 



( 2 . 6 ) 



where 

N = number of elements, 

= element excitation coefficient. 
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= kd X cos (0) +P , 



(2.7) 



and 

k = wave number, 

d = distance between elements (m) , 

0 = direction of main beam (degrees) , 

P = phase shift between elements (degrees) [Ref. 8] . 

The main beam can be scanned in any direction by 
adjusting the progressive phase shift between elements. 

M arrays of N elements can be arranged to form a 
planar array which will allow the beam to be scanned in any 
direction. Additionally, array elements can be excited 
non-uniformly to control the beam shape. If a coordinate 
system such as that of Figure 2.3 is chosen, it is shown in 
[Ref. 8] that the array factor in the X direction is given by 

M 

^ j ^j(i7?-l) (iccfj^in(0)cos{<|>) +Pjf) (2.8) 

m=l 

similarly, the array factor in the Y direction is given by 

N 

I +pp (2.9) 

75 = 1 

where 

= excitation coefficient for m elements and 
I-, = excitation coefficient for n elements. 

-‘•n 
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z 




Figure 2.3. Array Coordinate System 

The total array factor is the product of and AP^ 
and is given by 



AF = 



M 

E 

m=l 



^mie 



j (;n-l) (icdjfSin (0) cos (<()) +p^f) ' 



N 






jin-l) {kdySin{e)sin{iit) +Py) 



n=l 



( 2 . 10 ) 

2. Beam Synthesis 

The beam in a phased array antenna is shaped by non- 
uniform excitation of the array elements. This non-uniformity 
can be in amplitude, phase, or both. Using general synthesis 
[Ref. 9] and optimization [Refs. 10 and 11] techniques, beams 
of almost any shape and contour can be formed. For less 
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stringent beam shape criteria, a simpler synthesis method 
based on the Fourier transform method presented in [Ref. 8] 
has been derived. 

In the Fourier transform method, excitation 
coefficients are related to the Fourier transform of the array 
factor expressed as a function of y, which is itself a 
function of 0 and <{). Excitation coefficients are given by 






-N n i N 



( 2 . 11 ) 



for an odd number of elements, and 



= 



— f AF(\|;) e 2 
2 TC ~7c 



(2.12) 



-M < m -1 



and 



= 



^TAFiiir) e 

2Jt J -n 



-j 1 ] rif 



di|f 



1 < m < M 



(2.13) 



for an even number of elements [Ref. 8] . 

Equations 2.11 - 2.13 yield complex numbers which 

represent the amplitude and phase excitation of an element . 
The total excitation coefficient is obtained as the product of 
and J„. 

n m 
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Ill . ANTENNA DESIC3N 



A. MICROSTEaP PATCH DESIGN 
1 . S\jbstrate Material 

The first step in the design process is the selection 
of a dielectric siibstrate material . The dominant features of 
a microstrip array are controlled by substrate parameters such 
as thickness and permittivity more than by the particular 
element type [Ref. 12] . 

Because weight is a primary consideration in 
spacecraft design, a substrate must be chosen which not only 
has satisfactory dielectric properties, but which also has low 
density. There are a variety of substrate materials available 
with a wide range of relative dielectric constants and 
densities. Tables summarizing various materials and their 
relevant properties can be found in [Refs. 5 and 6] . 

Two representative substrate materials which were 
considered for this design are Duroid (RT 5870) , produced by 
Rogers Corporation of Chandler, AZ, and Kevlar (HI 4093) , 
produced by Arlon Electronic Substrates Division of Rancho 
Cucamonga , CA . 

Studies have shown [Ref. 13] that efficiency and 
bandwidth are nearly independent of patch shape ( i . e . , 
circular vs rectangular) and are determined mainly by 
substrate thickness and permittivity. 



18 



In order to achieve good resolution in beam shape and 
steering, the number of array elements must be maximized 
within available array area and element spacing constraints . 
This is accomplished by making the microstrip patches as small 
as possible while maintaining acceptable performance charac- 
teristics. For a given operating frequency, patch size is 
determined by relative dielectric constant and substrate 
thickness. Figure 3.1 shows the effect of relative dielectric 
on patch width based on equation (2.1) . Figure 3.2 shows the 
effect of dielectric thickness (h) on patch length for a given 
operating frequency and relative dielectric constant, based on 
equation (2.2) . 

Duroid substrates are available with ranging from 
2.1 to 10.0. Higher yields a smaller patch size at the 
expense of bandwidth and gain. Rogers Duroid RT 5870 has 
equal to 2.23 which allows higher efficiency and bandwidth 
performance but the density of RT 5870 Duroid is 2.2 gm/cm^. 
Arlon's HI 4093 Kevlar material on the other hand, has a 
relative dielectric constant of 3.9 and a density of 1.38 
gm/ cm^ . 

For a given operating frequency and dielectric 
thickness, Kevlar yields a smaller patch size than Duroid and 
although efficiency and bandwidth performance are reduced, the 
results are still within acceptable parameters. Therefore, 
the lower mass of Kevlar far outweighs its disadvantages. 
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Patch Width vs. Relative Dielectric 
(freq=1618 .25 MHz) 




Relative Dielectric ( ^r) 



Figure 3.1. Effect of on Width of a Microstrip Patch 
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Patch Length vs. Dielectric Height 
(freq=1618 .25 MHz) 




Figure 3.2. Effect of Height and on Length of a 

Microstrip Patch 
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Using a solid Duroid dielectric 0.3 cm thick, the mass 
of the substrate alone would be 36.7 kg (for an antenna area 
ofl.4mx4.0m).A Kevlar substrate with the same dimensions 
would yield an antenna mass of 23.3 kg. Appendix B contains 
detailed mass estimates for antenna components. 

2. Patch Dimensions 

a. Ihickness 

Because the antenna substrate and ground plane 
must form the array panels, they must be made thick enough to 
provide the required rigidity. The array panels must be rigid 
enough to withstand structural loads of launch and deployment 
as well as resist warpage due to thermal cycling. A siibstrate 
thickness of at least 0.3 cm was selected and assumed adequate 
to meet these criteria. The antenna mass could be reduced 
without significantly affecting electrical and mechanical 
performance by supporting a thin dielectric siibstrate on a 
honey-comb structure such as that described in [Refs. 3 and 
4] . The effect of raising the height of the patch above the 
ground plane is to increase bandwidth [Ref. 13] , while 
reducing dielectric thickness would have the opposite effect . 

b. Area 

The equations presented in Chapter II can be used 
as a starting point for determining patch dimensions which can 
later be optimized using a more sophisticated model . For an 
operating frequency of f^=1618.25 MHz and a relative 
dielectric constant e^=3.9, equations (2.1) -(2.4) yield a 
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patch length of 5.8 cm and a width of 7.2 cm (note that the 
width is defined as the long dimension of the rectangle) . The 
gain of a single patch is approximately 5 dB. Detailed cal- 
culations for length, width, efficiency, directivity, and gain 
calculations for a single patch are contained in Appendix C. 

A FORTRAN computer program which makes use of 
modal expansion theo 2 ry was used in an iterative process to 
optimize the values for patch length given a constant width. 
A listing of the program is contained in Appendix F. The 
final values for length and width are summarized in Table III- 
1 along with the operating parameters for a single microstrip 
patch. 



TABLE III-l. PARAMETERS FOR A SINGLE MICROSTRIP PATCH 



e =3.9, h=0.3 
(cm) 

f, =1618. 25 MHz 


L (cm) 


W (cm) 


G (dB) 


% BW 


Transmission Line 


4.6 


5.9 


5.02 


2.42 


Moda 1 Expans i on 


4.44 


5.9 


5.02 


2.62 



B. PLANAR ARRAY DESIGN 

Six array panels are needed to produce the six ellipti- 
cally shaped beams . The beams are symmetrical from outboard 
to inboard with beams 1 & 6 being identical, beams 2 & 5 
identical, and beams 3 & 4 identical. The total area avail- 
able for each array panel in the deployed configuration is 
0.9338 m^. 
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1 . Element Spacing 

As previously stated, the number of array elements 
needs to be maximized in order to achieve good beam 
resolution. At the same time, the array elements must be 
placed far enough apart to avoid mutual coupling which would 
degrade performance by increasing side lobe levels (SLL) . It 
has been shown [Ref. 14] that for an element spacing of d > 
'k/2, mutual coupling will not cause significant side lobe 
levels. Measurements made on L-band microstrip arrays have 
shown that mutual coupling levels of -25 dB or less can be 
achieved with proper element spacing [Refs. 14 and 15] . If 
the elements are separated by much more than one half 
wavelength, unwanted grating lobes are produced. An element 
spacing of about 0 . 55X has been shown to give good results 
[Ref. 3] . To avoid grating lobes and mutual coupling effects, 
the patch spacing for this design was chosen as 0.1017 m 
(0.55X.) . This spacing allows a 4 x 9 element array to be 
placed on each of the six antenna panels . 

2. Excitation Coefficients 

The Fourier transform method which was previously 
developed for a planar array was used to obtain element phase 
and amplitude excitation coefficients. Array coefficients for 
a single beam are shown in Table III-2. The resulting array 
factors are plotted in Figure 3.3 and Figure 3.4. The excita- 
tion coefficients are tabulated in Appendix D and array factor 
plots for each beam are presented in Appendix E. 
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TABLE III -2. BEAM 1 EXCITATION COEFFICIENTS 



Element 

Number 


1 


2 


3 


4 


1 


0.013 

Zl52.36° 


0.046 
L-121 .6° 


0.046 
L-121 .6° 


0.013 

LS2.2S° 


2 


0.019 

A-57.73° 


0.066 

^129.3° 


0.066 

L129.2° 


0.019 
L-Sl .12° 


3 


0.024 

Zl-153.8° 


0.083 

L2Q.1Q° 


0.083 

A28.18° 


0.024 
Zl-153 .8° 


4 


0.027 
yil03 .1° 


0.094 

L-1S.3 


0.094 

L-ie.9° 


0.027 
A103 .1° 


5 


0.028 

A0.0° 


0.098 

A.180.0° 


0.098 

A-180.0° 


0.028 

A0.0° 


6 


0.027 

A-103.1° 


0.094 
Lie. 9° 


0.094 

L16.9° 


0.027 
A-103 .1° 


7 


0.024 
/!l153 .8° 


0.083 

i-28.18° 


0.083 

A-28.18° 


0.024 

il53.8° 


8 


0.019 
LSI .12° 


0.066 

A-129.3° 


0.066 

L-129.2° 


0.019 
LSI .12° 


9 


0.013 

L-S2.26° 


0.046 
L121 .S° 


0.046 

A.127.6° 


0.013 

L-S2.26° 



3 . Antenna Gain 

The maximum array gain is the sum of individual 
element gains less losses associated with the array feed. The 
maximum gain is 20.6 dBi. This is well above the maximum 
required (6 dBi) for any antenna panel, however, a significant 
amotint of loss will occur in the feed structure and must be 
offset by the excess gain. 
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Figxire 3.3. Beam 1 Array Factor (<1)=90°) 




Figure 3.4. Beam 1 Array Factor (<{)=0°) 
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4. Feed Methods 



There are a variety of feed methods for microstrip 
patch arrays and for each method there exists many variations . 
In general, feed methods are grouped into three categories - 
probe, corporate, and triplate feeds - and are briefly 
discussed next . 

a. Probe Feed 

In the probe feed method, each patch is connected 
to the feed network via a coaxial -fed probe which extends 
through the ground plane and siibstrate. The probe is physi- 
cally attached to the patch by a solder joint. Impedance 
matching between the coax and array element can be accom- 
plished for thick substrates by designing a teardrop- shaped 
probe which cancels out probe inductance [Ref. 3] or by 
varying the diameter of the probe for thin substrates [Ref . 
16] . Obviously it would be impractical to construct an array 
where each element had its own coaxial feed cable. The weight 
and space requirements would be prohibitive for spacecraft 
applications . 

b. Corporate Feed 

A corporate feed offers a practical solution to 
the array feed problem by allowing one coaxial line to feed a 
network composed of microstrip transmission lines. A series 
of power dividers delivers power to and from the radiating 
elements. The microstrip transmission line provides equal 
path length to all elements maintaining phase between 
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elements. The two major disadvantages of the corporate feed 
are that there is not enough room on a planar array substrate 
for radiating elements and the feed network, and high line 
loss occurs due to long path lengths . The second disadvantage 
can be overcome by integrating active elements into the 
substrate [Refs. 17 and 18] or by breaking the feed network 
into subsections, each of which is fed by a low loss coaxial 
cable [Ref. 19] . 

c . Triplate Feed 

In a triplate feed (also called a hybrid feed) , 
the feed network is etched on its own substrate layer and 
either laminated to the underside of the radiating element ' s 
substrate or attached to the backside of the ground plane. 
The former method requires the radiation of the feed network 
to be accounted for in the antenna design while the latter 
method eliminates feed radiation effects and divides the 
design problem into two parts: the array itself and the feed 
network . 

(1) Feed Network in Front of Ground Plane. A 
triplate feed located in front of the grotind plane and 
directly under the radiating element's substrate simplifies 
the array construction at the expense of increased design 
complexity. The metal pins used to attach the array element 
to the feed network are eliminated through electromagnetic 
coupling. Impedance matching is controlled by the amount of 
overlap between the microstrip patch and the feed network. 
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Radiation effects can be largely reduced by locating the feed 
network as close to the grotind plane as possible. [Ref. 13] 
(2) Triplate Behind Ground Plane. If the 
corporate feed network is located behind the ground plane, 
radiation effects from the feed network are eliminated at the 
expense of a more complicated and expensive physical construc- 
tion involving plated through holes and connecting pins. 

5 . Feed Design 

Although antenna feed design is discussed in [Ref . 
20], a complete design is outside the scope of this thesis. 
A brief discussion of a possible feed method is included to 
prove the feasibility of providing the required array element 
excitations . 

A triplate corporate feed located behind the antenna 
ground plane could be used. This feed method, based on those 
presented in [Refs. 6 and 19] , would simplify antenna perform- 
ance analysis by eliminating feed radiation effects. Since 
the antenna is a statically scanned array, the required phase 
shifts and impedance matching could be built into the 
microstrip feed lines using microstrip transmission line 
techniques. Figure 3.5 shows one possible feed layout. 

C. SUMMARY 

An exploded view of one antenna panel is illustrated in 
Figure 3.6. 
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Figure 3.6. Antenna Panel Assembly 
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IV. PERFORMANCE ANALYSIS 



A. ANTENNA MODEL 

Chapter II described the wire grid, transmission line, and 
modal expansion models which are used to analyze and predict 
the performance of microstrip antennas. The wire grid model 
is unsuitable for use on any but the simplest of problems due 
to the computational resource required. Because the trans- 
mission line model and the modal expansion model do not 
account for surface waves in the dielectric substrate and do 
not account for mutual coupling, they do not provide the most 
accurate solution available [Ref. 16] . However, it will be 
shown that for arrays with proper element spacing, and all but 
extreme scan angles, the modal expansion modal provides 
adequate results with relative ease of computation. The areas 
where the modal expansion model fails will be discussed and 
shown to be insignificant for this analysis. 

1. Infinite vs. Finite Arrays 

The analysis of an infinite array yields adequate 
results for the inner array elements and accounts for a change 
in impedance with varying scan angle. Infinite array analysis 
ignores edge effects due to surface waves in the dielectric. 
The analysis of a finite array includes these edge effects but 
at the expense of more difficult calculations. It has been 
shown that surface wave power increases with dielectric 
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thickness and that for a thin dielectric (h < 0.03A,) , surface 
wave power is negligible [Ref. 14]. It has also been shown 
that surface wave power diminishes with increasing array size 
and that for array sizes of 7 x 7 or greater an infinite 
approximation is acceptable [Ref. 21] . 

2 . Edge Effects 

Diffractions cause by the finite edge of the ground 
plane play a significant part in the radiation pattern at wide 
scan angles and in the backlobe region. The Geometric Theory 
of Diffraction (GTD) can be combined with other analysis 
methods to provide an accurate antenna performance model in 
these regions. Diffraction edge effects can be ignored for 
scan angles less than about 70° [Ref. 22] . 

3 . Input Iirpedance 

Although the modal expansion model enables calculation 
of input impedance, patch currents near the feed point are not 
modeled. This results in inaccurate impedance. It has been 
shown that for substrate thicknesses of less than about 0.02X., 
the idealized feed model provides adequate results [Ref. 23] . 

B. SOFTWARE MODELS 

1. ^thod of Moments 

The method of moments solves for \mknown currents on 
the surface of the microstrip patch. If a dielectric is 
present, the dielectric Green's function is used to solve the 
electric field integral equation. This method automatically 
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accounts for mutual coupling and provides accurate results for 
antenna performance. The method of moments depends on the 
accurate calculation of the elements in the impedance matrix 
for its precision [Ref. 24] and requires a large amount of 
computer resource . Several methods have been developed to 
calculate method of moment solutions more efficiently, [Refs. 
25 and 26] , but comparisons with results from other models do 
not show significant differences [Ref. 24] . 

2. I^lbdal Expansion Model 

The modal expansion model was discussed in Chapter II. 
The limitations of this model have been shown in previous 
discussions to be insignificant for this application. Because 
of its simplicity, the modal expansion model was used to 
predict the performance of this antenna. The software program 
used for this analysis was written in FORTRAN and obtained 
from the Jet Propulsion Laboratory in Pasadena, CA. A source 
code listing for this program is contained in Appendix F. 

C. MODELING SCENARIO 

Using the excitation coefficients previously obtained, 
antenna patterns for each of the six beams were generated for 
scan angles of 0=-9O® to +90®, <])=0 and 0=-9O® to +90®, (])=90. 
The results were plotted in order to asses whether or not each 
beam met the design specification in terms of scan angle and 
beam width. A contour plot of the approximate grovind 
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footprint of each beam was also generated in order to assess 
antenna ground coverage . 

D. RESULTS 

1 . Radiation Patterns 

The radiation pattern obtained for beam 1 is plotted 
in Figure 4.1. Beam patterns for the remaining beams are 
contained in Appendix G. 

The model results demonstrate that both beam shape and 
beam orientation are attainable for beam 1 . Similar results 
were obtained for beams 2 through 6 with one notable 
exception. The specification called for beams 2 and 5 to be 
scanned 30 degrees from center. When the beams were steered 
to exactly -30 degrees and 30 degrees respectively, the result 
resembled a broad side array beam pattern. Examination of the 
antenna's normal beam pattern reveals a null which prevents 
this antenna from being scanned to exactly 30 degrees. By 
adjusting the scan angle to -27 degrees for beam 2 and 27 
degrees for beam 5, acceptable beam patterns were obtained. 
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Figure 4.1. Beam 1 Radiation Pattern ((|)=90°) 



2 . Pattern Separation 

Recall that the beams are required to operate in 
transmit /receive pairs such that the satellite uplink and 
downlink can be operated simultaneously at the same frequency. 
This requires adequate pattern separation to avoid mutual 
interference at the earth ' s surface . Figure 4 . 2 demonstrates 
that this design achieves adequate pattern separation for 
beams 1 and 4 . Similar results were obtained for beam pair 2 
and 5 and pair 3 and 6 and are presented in Appendix G. 
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Figure 4.2. Simultaneous Beam Patterns 



3 . Antenna Footprint 

In order to clearly demonstrate antenna coverage, the 
relative gain for each beam was computed over a 5 degree by 5 
degree grid centered about the satellite's nadir. The results 
for beam 2 are presented in Figure 4.3. The beam pattern is 
shown in terms of scan angle from nadir in the in- track and 
cross-track directions. Similar results were obtained for the 
remaining beams and are contained in Appendix G. 
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Relative dB Contours 
(Beam scanned 27 degrees from nadir) 




Left Right 

Cross Track 



Figure 4.3. Beam 2 Relative dB Gain Footprint 
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V. ALTERNATIVE DESIC2^ 



A. ANTENNA CONFIGURATION 

An alternative design to that presented in the previous 
chapters is one which takes advantage of the fact that the L- 
band communications signal uses time division multiplexing 
(TDM) . The use of TDM, as implemented in the payload for this 
spacecraft, means that each of the six antenna beams has a 60 
ms duty cycle. Since adjacent pairs of antennas are not 
operating at the same time, the antenna could be constructed 
such that there are two beams per antenna panel thus reducing 
the number of panels required. 

An integrated multiple beam microstrip array similar to 
that described in [Ref. 27] could be constructed. Two feed 
lines would be constructed on a lower substrate which would be 
electromagnet ically coupled to the patch array on the 
substrate above. Each feed line would be energized at the 
appropriate time corresponding to the on-time for the 
particular beam. This configuration would simplify antenna 
construction by eliminating the need for through-hole plating 
and connecting pins between the feed network and the patch 
elements [Ref. 28] . The deployment mechanism for this 
configuration would also be much simpler as only two hinged 
panels would have to unfold as opposed to four. There would 
also be a reduction in weight while performance would actually 
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be enhanced due to reduced losses in the smaller feed network 



required. 

Combining two antenna panels into one and using that 
single panel for two beams would almost double the number of 
patch elements available to each beam. In addition to 
increasing gain, this would allow better resolution in both 
beam shape and orientation. Increasing the number of array 
elements would also eliminate unwanted nulls which adversely 
affect scan angle. 

B. ANTENNA ORIENTATION 

The antenna in this study is oriented such that individual 
panels are perpendicular to the spacecraft's velocity vector. 
This orientation was chosen for structural and thermal control 
reasons during the early stages of the preliminary design and 
is not optimum. If the design described above were adopted, 
the antenna panels could be oriented parallel to the velocity 
vector allowing cross track beam steering to be accomplished 
by tilting two outside panels towards the center. All 
available patch elements could then be dedicated to shaping 
the beam which would result in reduced beam overlap. 
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VI. SUMMARY 



A. ANALYSIS AND DESIGN 

Two different siibstrate materials were considered for the 
design of a multiple beam, phased array, microstrip antenna. 
Design curves were presented to illustrate tradeoffs between 
weight and performance in the substrate material. This 
information was used to choose an appropriate stibstrate 
material which would provide adequate performance at a minimum 
cost in mass . 

Closed form formulae were used to arrive at a preliminary 
width and length for microstrip patch elements. A computer 
model was then used to optimize the width and length 
dimensions . Array theory was applied to derive patch 
excitation coefficients. 

B. PERFORMANCE EVALUATION 

Patch excitation coefficients were used as input to a 
computer model, which employs modal expansion theory, and 
antenna radiation patterns were generated. 

While the Method of Moments would have yielded more 
accurate results, it would have been at the expense of 
increased computer time. The Modal Expansion method used 
provided adequate results with a reasonable amount of effort. 



40 



The results demonstrate that the proposed design would not 
meet all the design goals. While the desired beam shapes were 
obtained, the ground footprint of the outside beams (beam 1 
and 6) show that most of the energy would be concentrated 
outside the satellite's field of view. 

C. CX)NOiUSIONS 

Analysis showed that this antenna could not be constructed 
within the allocated mass budget. In fact, the mass of the 
substrate material alone would be over two times the allocated 
mass budget . 

The use of rectangular microstrip patches is not an 
efficient use of the available array area. An octagonally 
shaped patch would allow more elements to be placed on a 
panel. This would allow tighter control over beam shape and 
orientation . 

The alternative design discussed in Chapter IV would 
address the short -comings in the original design and allow the 
ground footprint to be entirely within the satellite ' s field 
of view. However, operating two beams per antenna panel could 
require active phase shifters, increasing cost as well as 
reliability risk. 

The use of fiber optic cables instead of coaxial transmis- 
sion line and etched copper feed networks would reduce both 
weight and power loss. 
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since the methods for determining patch dimensions 
discussed in the literature all give slightly different 
results, this antenna would have to be tested and trimmed 
after construction in order to ensure the operating parameters 
were met . 

D. FUTURE TEIESIS OPPORTUNITIES 

A scale model could be built and tested by increasing the 
operating frequency to yield a smaller wavelength. The 
coupled feed model discussed in the alternative design could 
also be built and tested. Various orientations, beam 
patterns, and scan angle could be explored yielding useful 
data for future phased array design projects. 
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APPENDIX A 



LINK BUDGETT SUMMARY 

TABLE A-1. LINK BUDGET [REF. 1] 





User to 
Satellite 






Satel 1 i te 
to User 


t 

t 


Frequency 


1625.0 


MHz 


Frequency 


1625.0 


MHz 1 


RF Power 


6 


Watts 


RF Power 


19 


Watts 




7.8 


dBw 




12.8 


dBw 


Power Loss 


-1.0 


dB 


Power Loss 


-2.2 


db 


Antenna Gain 


3.0 


DBi 


S/C Ant. Gain 
( isoflux) 


4. 0-6.0 


dBi 


EIRP 


9.8 


dBw 


EIRP 


14.6 


dBw 














Satel 1 ite 
Altitude- 
750 NM 


1389.0 


km 


Satellite 

Altitude 


1389.0 


km 


Elevation Angle 


90.0 


degrees 


Elevation 

Angle 


90.0 


degrees 


Range 


U25.0 


km 


Range 


1389.0 


km 


Free Space Loss 


-159.5 


dB 


Free Space 
Loss 


-159.5 


dB 


RX Signal 
Strength 


-149.8 




RX Signal 
Strength 


-159.5 




Polarization 

Loss 


-0.5 


dB 


Polarization 

Loss 


-0.5 


dB 


Tracking Loss 


0.0 


dB 


Tracking Loss 


0.0 


dB 


S/C Ant. Gain 
(isoflux) 


4. 0-6.0 


dBi 


Antenna Gain 


3.0 


dBi 


RX Line Losses 


-1.0 


dB 


RX Line Loss 


-0.5 


db 














L-Band Amplifier 
Eff. 


35 


% 


L-Band 

Amplifier Eff. 


35 


% 


Ej,/NoRequi red 


3.5 


dB 


Eij/Nq Required 


3.5 


dB 


Total BU 


16.5 


MHz 


Total BU 


16.5 


MHz 
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APPENDIX B 



I . MASS ESTIMATES 

A. Kevlar Sijbtrate 

Density = 0.05 Ib/in^ Density = 1.4 gm/cm^ 

The total siibstrate volume is given by: 

1 = 1.4 m w = 0.667 m h = 0.003 m 
Voliome = lxwxhx2x6 
Volume = 3.4 x 10^ cm^ 

The total substrate mass is given by: 

Mass,(gy^g^ = Volume x Density 

Mass,(gy^g^ = 46.5 Kg Mass,(g^^g^ = 102.6 lb 

B. Duroid Substrate 

Specific_gravity = 2.2 Density = 2.2 gm/cm^ 

Total substrate mass: 

MasSp^j^gj^ = Volume x Density 

= 74 Kg MasSp^^g,.^ = 163 lb 
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C. Copper 

1. Microstrip Patches 

Densityj.y = 8.96 gm/cm* 

Total number of elements =9x4x6= 216 
Element width = 5.9 cm Element length = 4.6 cm 

ThicknesSj-^j = 0.0127 cm 

Total volume of array elements = length x width x thickness^.^ 
X tot num elements 

Total volume = 74.5 cm^ 

Total MasSg^g^g^^g = Density x volume 

Total Mass ^^^^„^3 = 0.7 Kg 

2. Microstrip Feed Lines 

lengthfggj = 1.3 m width^^^^ = 0.006 m 

Volume^ggj = length^ggj x width^^^^ x thickness^.^ 

= I^ensityj,^ X Volume^^^^ 

MasSfggd = 8.9 X 10*^ Kg 

3. Total Copper Mass 

The total copper mass is given by: 

Massc^ = Mass^i^^^,^ + Mass^^^d 
MasSj,^j = 0.701 Kg 

D. Muminum (^x>tind Planes 

Since a shielded triplate feed is used, two aluminum 
ground planes are needed. 

Density^^ = 2.7 gm/cm^ thickness^^ = 1.0 x 10’^ 

Volume^^ = 1 X w X thickness^^ x 4 
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Mass^^ = Volume^ ^ x Density^^ x 2 
Mass^j = 1.35 X 10'^ Kg 



II. SUMMARY 

The Duroid substrate is approximately 1.5 times heavier 
than the Kevlar for a given volume. Therefore, the Kevlar 
material should be used for the substrate. The total mass for 
the antenna is estimated as follows: 

Mass = Mass,^gy^g^ + Mass^.^ + Mass^^ 

A mass margin of 0.01 is chosen and the total antenna mass 
becomes 



Mass 



antenna 



Mass X (1 + Mass, 



'margin 



Mass 



antenna 



52 Kg Mass 



antenna 



114.5 lb 



46 



APPENDIX C 



I . PATCH DESIGN 

Using the transmission line model, the patch parameters 
are calculated as follows : 

Given: f = 1618.25 MHz = 3.9 h=0.3 cm 

1 = 3(10)Vf k = 2 tiA 

A. = 18.5 cm k = 33.9 






c 

2fr 



2 



2 

2 



(cm) 



W = 5.9 cm 



The effective dielectric constant is given by: 



Cr+l 



e.-l, 



1 + 



12h\- 
Pi I 



e„ = 3.6 

e 



The patch length is then calculated as follows : 

(e^+0.3) (-^+0.264) 

AJ = 0.412h A1 = 0.1 cm 

(e^-0.258) (-^+0.8) 

L = — - 2AJ (cm) 

L = 4.6 cm 
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II. GAIN 



In order to calculate gain, the directivity must be 
calculated according to the following: 



^0 = 






|gin( )2tan(0)2sin(0)de 



Losses, expressed as resistance, must also be found according 

to the following relations : 

Quality factor: 

_ 3(10)^/E; 

^ 4 fh 



Radiation resistance: 
Copper loss : 



= 120A/W 

= 2 .7 (10) 






(-^) {%) 

10^ ^ 



Dielectric loss: 



= 



_ 81 (10) h\ 



L W 



The total resistance is then R^ = R^, + R^, + R^j ^ = 379.8 
With a of VSWR = 2, Bandwidth (BW) and efficiency (t|) are 
found as follows: 



BW = X 100 

{Q sJvmR) 



Directivity and Gain are computed as : 
D = 10 X log(Do) 

G = 10 X log(T| X Dq) 



n = 



Rn 



D = 5.06 dB 
G = 5.0 dB 
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APPENDIX D 



EXCITATION COEFFICIENTS 



TABLE D-1. BEAM 1 EXCITATION COEFFICIENTS 



Element 

Number 


1 


2 


3 


4 


1 


0.013 

/i.52.36° 


0.046 

4.-127.6° 


0.046 

4.-127.6° 


0.013 

4.52.36° 


2 


0.019 

A-57.73° 


0.066 

4.129.3° 


0.066 

4.129.3° 


0.019 

4l-57.73° 


3 


0.024 

/i-153.8° 


0.083 

4.28.18° 


0.083 

4.28.18° 


0.024 

4.-153.8° 


1 


0.027 
^103 .l" 


0.094 

A-76.9° 


0.094 
4.-76. 9° 


0.027 
4.103 .1° 


5 


0.028 

4.0.0° 


0.098 

4l180.0° 


0.098 

4.-180.0° 


0.028 

4.0.0° 


6 


0.027 
4.-103 .1° 


0.094 

4.76.9° 


0.094 

4.76.9° 


0.027 
4.-103 .1° 


7 


0.024 

4.153.8° 


0.083 

4.-28.18° 


0.083 

4.-28.18° 


0.024 
4.153 .8° 


8 


0.019 

4.57.73° 


0.066 

4l-129.3° 


0.066 

4.-129.3° 


0.019 

4l57.73° 


9 


0.013 

4.-52.36° 


0.046 

4.127.6° 


0.046 

4.127.6° 


0.013 

4.-52.36° 
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TABLE D-2. BEAM 2 EXCITATION COEFFICIENTS 



Element 

Number 


1 


2 


3 


4 


4 


0.003 
Lies. 2° 


0.010 

4.-10.68° 


0.010 

4.-10.68° 


0.003 

4.169.3° 


2 


0.011 

A.-55.0° 


0.039 

4.127.0° 


0.039 

4.127.0° 


0.011 

4l-55.0° 


3 


0.029 

/i-95.3° 


0.099 

4l84.6° 


0.099 

4.84.6° 


0.029 

4l-95.3° 


4 


0.043 

L-137.7° 


0.148 

4.42.3° 


0.148 

4.42.3° 


0.043 

4.-137.7° 


5 


0.048 

Zl-180.0° 


0.167 

4.0.0° 


0.167 

4.0.0° 


0.048 

4.180.0° 


6 


0.043 

A.137.7° 


0.148 

4l-42.3° 


0.148 
4.-42. 3° 


0.043 

4.137.7° 


7 


0.029 

A.95.3° 


0.081 
4.-84. 6° 


0.099 
4.-84. 6° 


0.029 

4.95.3° 


8 


0 . Oil 
4.53 .0° 


0.039 

4.-127.0° 


0.039 

4.-127.0° 


0.011 
4.53 .0° 


9 


0.003 

4.-169.3° 


0.010 

4.10.68° 


0.010 

4l10.68° 


0.003 

4l-169.3° 
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TABLE D-3. BEAM 3 EXCITATION COEFFICIENTS 



Element 

Number 


1 


2 


3 


4 


4 


0.013 
LI .1° 


0.046 

A-172.3° 


0.046 
4.-172 .3° 


0.013 

4.7.7° 


2 


0.002 

Lb.Q° 


0.007 

4.-174.2° 


0.007 
4.-174 .2° 


0.002 

4.5.8° 


3 


0.024 


0.083 
4.3 .8° 


0.083 
4.3 .8° 


0.024 
4.-176 .1° 


4 


0.051 

4.-178.1° 


0.176 

4.1.9° 


0.176 

4.1.9° 


0.051 

4.-178.1° 


5 


0.062 

4.-180.0° 


0.215 

4.0.0° 


0.215 

4.0.0° 


0.062 

4.180.0° 


6 


0.051 

4.178.1° 


0.176 
4.-1. 9° 


0.176 
4.-1. 9° 


0.051 

4.178.1° 


7 


0.024 

4l176.1° 


0.083 
4.-3. 8° 


0.083 
4.-3 .8° 


0.024 

4.176.1° 


8 


0.002 
4.-5. 8° 


0.007 

4.174.2° 


0.007 

4l174.2° 


0.002 
4.-5. 8° 


9 


0 . 013 
L-1 .1° 


0.046 
4.172 .3° 


0.046 
4.172 .3° 


0.013 
4.-7. 7° 
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TABLE D-4. BEAM 4 EXCITATION COEFFICIENTS 



Element 

Number 


1 


2 


3 


4 


4 


0.013 
L-1 .1° 


0.046 

4.172.3° 


0.046 

4.172.3° 


0.013 
4.-7. 7° 


2 


0.002 

A.-5.8° 


0.007 

4l174.2° 


0.007 

4.174.2° 


0.002 
4.-5. 8° 


3 


0.024 

A.176.1° 


0.083 
4l-3 .8° 


0.083 
4l-3 .8° 


0.024 

4.176.1° 


4 


0.051 

L.11Q.1 


0.176 

4l-1.9° 


0.176 
4.-1. 9° 


0.051 

4.178.1° 


5 


0.062 

^-180.0° 


0.215 

4.0.0° 


0.215 

4.0.0° 


0.062 

4.180.0° 


4 


0.051 

4.-178.1° 


0.176 

4.1.9° 


0.176 

4.1.9° 


0.051 

4.-178.1° 


7 


0.024 

4.-176.1° 


0.083 
4.3 .8° 


0.083 

4l3.8° 


0.024 

4.-176.1° 


8 


0.002 

4.5.8° 


0.007 

4.-174.2° 


0.007 
4.-174 .2° 


0.002 

4.5.8° 


9 


0.013 
LI .1° 


0.046 

4.-172.3° 


0.046 
4.-172 .3° 


0.013 

4.7.7° 



52 



TABLE D-5. BEAM 5 EXCITATION COEFFICIENTS 



Element 

Number 


1 


2 


3 


4 


1 


0.003 

A.-169.3° 


0.010 

A.10.68" 


0.010 

4l10.68° 


0.003 

4.-169.3° 


2 


0.011 

/<l55.0” 


0.039 

^-127.0° 


0.039 

4.-127.0° 


0.011 

4.55.0° 


3 


0.029 

L95.3° 


0.099 

L-84.6° 


0.099 

A-84.6° 


0.029 

4.95.3° 


1 


0.043 

A.137.7° 


0.148 
4.-42. 3° 


0.148 

4l-42.3° 


0.043 

4.137.7° 


5 


0.048 

Zl-180.0° 


0.167 

4.0.0“ 


0.167 

4.0.0° 


0.048 

4.180.0° 


6 


0.043 

A.-137.7° 


0.148 

4l42.3° 


0.148 

4.42.3° 


0.043 

4.-137.7° 


7 


0.029 

A-95.3° 


0.081 
4l84 . 6“ 


0.099 

4.84.6° 


0.029 
4.-95. 3° 


8 


0.011 

^-53.0° 


0.039 

4.127.0° 


0.039 

4.127.0° 


0.011 
4.-53. 0° 


9 


0.003 

Zl169.3° 


0.010 

4.-10.68° 


0.010 

4.-10.68° 


0.003 

4.169.3° 
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TABLE D-6. BEAM 6 EXCITATION COEFFICIENTS 



Element 

Number 


1 


2 


3 


4 


1 


0.013 
A.-52 .36° 


0.046 

^127.6° 


0.046 

4.127.6° 


0.013 

4.-52.36° 


2 


0.019 

4.57.73° 


0.066 

4l-129.3° 


0.066 

4.-129.3° 


0.019 

4.57.73° 


3 


0.024 

4.153.8° 


0.083 

4l-28.18° 


0.083 

4.-28.18° 


0.024 

4l153.8° 


3 


0.027 
4.-103 .1° 


0.094 
Lie. 9 


0.094 

4.76.9° 


0.027 

4.-103.1° 


5 


0.028 

4l0.0° 


0.098 

4.180.0° 


0.098 

4l-180.0° 


0.028 

4.0.0° 


6 


0.027 

4.103.1° 


0.094 
4.-76. 9° 


0.094 
4.-76. 9° 


0.027 
4.103 .1° 


7 


0.024 

4.-153.8° 


0.083 

4.28.18° 


0.083 

4.28.18° 


0.024 

4l-153.8° 


8 


0.019 

4.-57.73° 


0.066 

4.129.3° 


0.066 

4.129.3° 


0.019 

4.-57.73° 


9 


0.013 
4.52 .36° 


0.046 

4.-127.6° 


0.046 

4.-127.6° 


0.013 

4l52.36° 
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APPENDIX E 



ARRAY FACTOR PLOTS 




Figiore E-1. Beam 1 Array Factor ((|)=90°) 




Figure E-2. Beam 1 Array Factor (<j)=0°) 
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Figiire E-4. Beam 2 Array Factor (<|)=0”) 
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Figure E-5. Beam 3 Array Factor (4>=90) 




Figure E-6. Beam 3 Airray Factor (4>=0‘’) 
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I 




Figure E-7. Beam 4 Array Factor ({j)=90°) 
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Figure E-9. Beam 5 Array Factor ((j)=90°) 




Figure E-10. Beam 5 Array Factor (<J>=0‘’) 
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Figure E-11. Beam 6 Array Factor (4)=90‘’) 




Figure E-12. Beam 6 Array Factor (4)=0'’) 
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APPENDIX F 



I. DESIGN. FOR 



C. . .MICROSTRIP ANTENNA DESIGN CHARACTERISTICS CALCULATED BY 
MODAL 

C... THEORY, ALL INPUT LENGTH ARE IN CENTIMETER, BY DR. HUANG 
COMPLEX CJ, YWY,AY,DEL(5) , Cl, C2 , D4 , SKY, FC, WC 
REAL A,B,EPS 
CJ= (0. , 1. ) 

PI=3 . 14159265 

TPI=2.*PI 

DPR=180./PI 

TOL=.0001 

WRITE(6, 11) 

11 FORMAT (' CENTER FREQ=?') 

READ (5, 99) FREQ 
WRITE (6, 1) 

1 FORMAT (' WIDTH OF PATCH=? ' ) 

READ(5, 99)A 

99 FORMAT (F5. 8) 

WRITE(6,2) 

2 FORMAT (' RESONANT LENGTH OF PATCH=? ' ) 

READ(5, 99)B 

WRITE (6, 3) 

3 FORMAT (' DIELECTRIC CONSTANT=? ' ) 

READ (5, 99) EPS 

WRITE (6, 4) 

4 FORMAT ( ' DIELECTRIC THICKNESS=? ' ) 

READ(5, 99)T 

WRITE (6, 5) 

5 FORMAT ( ' PROBE LOCATION=? ' ) 

READ (5, 99) YO 

WRITE (6, 6) 

6 FORMAT ( ' LOSS TANGENT=? ' ) 

READ(5, 99)TLOS 

40 FF=1./SQRT(1.+10.*T/A) 

EPE= (EPS+1 . ) /2 . + (EPS-1 . ) /2 . *FF 

FF= (A/T+0.262) /(A/T+0.813) * (EPE+0 .3) / (EPE-0 .258) *0 .412 
WAMDA=2 . 1*B*SQRT(EPS) 

GW=0 . 00836*A/WAMDA 

BW=0 . 01668*FF*A*EPE/WAMDA 

FY=0. 7747+0. 5977* (A/B-1. ) -0.1638* (A/B-1. ) **2 . 

YWY=CMPLX (GW, BW) 
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AY=CJ*TPI*376 . 7/WAMDA*T/A*YWY*FY*B 
DEL(l) = (0. , 0. ) 

DO 20 1=1,4 

Cl=2 . *AY* (CMPLX (PI , 0 . ) -DEL ( I ) ) 

C2=AY*AY+2.*DEL(I) *PI-DEL(I) *DEL(I) -CMPLX (PI*PI , 0 . ) 
DEL(I+1) =C1/C2- (DEL (I) ) **3 ./3 . 

20 CONTINUE 
D4=DEL(5) 

SKY= ( CMPLX ( PI, 0.) -D4)/B 
FC=CMPLX (1 . , 0 . ) -CJ*TLOS 
WC=3 . ElO/SQRT (EPS*FC) *SKY 
FR=REAL(WC) /TPI 
QR=REAL (WC) /2 . /AIMAG (WC) 

FF=COS (PI*YO/B) 

CC=EPE*8 . 854*1 . E-14*A*B/2 . /T/FF/FF 
RAD=QR/REAL (WC) /CC 

RCU=0 . 00027*SQRT (FR/l .E+9) *B/A*QR*QR 

RDI=30 . *TLOS*T*WAMDA/EPS/A/B*QR*QR 

RIN=RAD+RCU+RDI 

Q=RIN*REAL (WC) *CC 

DFR=FR/Q 

EFF=RAD/RIN 

FR= FR/l. E+9 

DFREQ=FREQ-FR 

IF( (ABS (DFREQ) ) .LE.TOL)GO TO 42 
WRITE ( 6 , 4 1 ) FR , DFREQ 
IF(FR.LT.FREQ) THEN 
B=B- .0001 
ELSE 

B=B+.0001 
END IF 
GOTO 40 



41 FORMAT ( ' FREQUENCY= ' , FIO . 5 , 3X, ' DFREQ= ' , FIO . 5 ) 



42 WRITE (6, 31) FR 

31 FORMAT (' RESONANT FREQUENCY IN GIGAHERTZ= ' , FIO . 5) 
WRITE(6,32)RIN 

32 FORMAT ( ' INPUT RESISTANCE= ' , FIO . 5 ) 

DFR=DFR/l.E+6/2 . 

WRITE(6,33)DFR 

33 FORMAT ( ' BANDWIDTH IN MEGAHERTZ= ' , FIO . 5 ) 
WRITE(6,34)EFF 

34 FORMAT ( ' RADIATION EFFICIENCY= ' , FIO . 5 ) 

PERCBW=1/ (Q*SQRT (2.0)) 

WRITE(6,37)PERCBW 
FORMAT ( ' % BW= ' , FIO . 5 ) 



37 
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WRITE ( 6 , 3 8 ) RCU , RDI , RAD , Q 

38 FORMAT (' RCU= ' ,F10.5,2X, 'RDI=' ,F10.5,2X, ' RAD= ' ,F10.5, 
2X, 'Q=' ,F10.5) 

WRITE (6, 39) A, B 

39 FORMAT ( 'WIDTH=' ,F10.5,5X, ' LENGTH= ' ,F10.5) 

STOP 

END 



II. MOiCARY.FOR 



C. . .LINEARLY POLARIZED PLANAR ARRAY FORMED BY LINEARLY 
C. . .POLARIZED MICROSTRIP ELEMENTS 

C...WITH CROSS-POL, MODAL EXPANSION MULTIMODES ARE USED. 
C...NO GTD DIFFRACTION, ALL DIMENSIONS IN WAVELENGTHS 
C. . .PROGRAMMED BY DR. J. HUANG AT JET PROPULSION LAB. 

DIMENSION DBR(362) ,DBL(182) ,XAX(362) ,XX(2) ,YY(2) , 
HANG (40) 

DIMENSION PHASR(182) ,PHASL(182) ,XE(40) , YE(40), 
1AMP(40), PHAS(40) 

COMPLEX CJ, CJJ, EX, EY, EZ , XO , YO , ZO , ETH, EPH, ER, EL 
COMPLEX FX,FY,FZ,X1,Y1,Z1,X2,Y2,Z2,X3,Y3,Z3 
COMMON/DDl/EPS , T , A, B , 12 , IP , XS , YS , ZS , CJJ 
COMMON/DD2/CJ, PI , TPI , DPR 
COMMON/DD4/THPR, PHPR, RR 

C. . .MODIFIED TO READ DATA FROM A GEOMETRY FILE 
C...BY BILL BARFIELD, NPGS 
OPEN (33, FILE='DATA') 

OPEN (35, FILE= ' GEOMETRY ' ) 

CJ=(0. ,1.) 

CJJ= (l.E-15, l.E-15) 

PI=3 .14159265 
TPI=2.*PI 
DPR=180 . /PI 
WRITE (6,1) 

1 FORMAT (' DIELECTRIC CONSTANT AND THICKNESS=? ' ) 

C READ (5, 99) EPS, T 

READ (33, *) EPS, T 
WRITE (6, 99) EPS,T 

99 FORMAT ( F8.6,F8.6) 

98 FORMAT (18) 

WRITE (6, 2) 

2 FORMAT (' MICROSTRIP WIDTH AND LENGTH=? ' ) 
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C READ(5, 99) A,B 

READ(33, *)A,B 
WRITE(6, 99)A,B 

. . .A,B ARE THE PHYSICAL WIDTHS OF THE METALLIC PATCH 

A ALONE X-AXIS, B ALONE Y-AXIS, X-AXIS IS HORIZONTAL 

WRITE (6, 3) 

3 FORMAT ( ' FIELD POINT DISTANCE AND PATTERN CUT 
/ANGLE=? ' ) 

C READ(5,99)RF,PHD 

READ(33, *)RF, PHD 
WRITE (6, 99) RF, PHD 
WRITE(6,*)RF,PHD 

WRITE (6, 5) 

5 FORMAT (' FEED PROBE DIAMETER & DISTANCE FROM EDGE=? ' ) 

C READ(5,99)D,YP 

READ(33, *)D, YP 
WRITE(6, 99)D,YP 

WRITE (6, 4) 

4 FORMAT (' NUMBER OF ELEMENTS=? ' ) 

C READ(5,98)N 

READ(35, *)N 
WRITE ( 6 , * ) N 

WRITE(6, 6) I 

6 FORMAT (' NUMBER' ,12, 2X, 'X,Y LOCATION, AMPLITUDE, PHASE & 
/ORIENTATION=? ' ) 

DO 105 1=1, N 

C READ(5, *)XW, YW,AMG,DEG, lANG(I) 

READ(35,*)XW, YW,AMG,DEG, lANG(I) 

WRITE(6, *)XW, YW,AMG,DEG, lANG(I) 

XE(I)=XW 
YE(I)=YW 
AMP ( I ) =AMG 
PHAS ( I ) =DEG/3 6 0 . *TPI 
105 CONTINUE 
D=5.*D 

SK2=EPS*TPI*TPI 
PIB=(PI/B) **2 
PIA=(PI/A) **2 
Fl=l . /A/B/SK2 
Fl=l.E-7 
CYP=COS(PI*YP/B) 
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F2=2 . /A/B/ (SK2-PIB) *CYP 
SD=SIN(PI*D/A) / (PI*D/A) 

F3=4 . /A/B/ (SK2-4 . *PIA-PIB) *CYP*SD 
CYP=COS (TPI*YP/B) 

F4=4 . /A/B/ (SK2-4 . *PIA-4 . *PIB) *CYP*SD 

F5=2 . /A/B/ (SK2-4 . *PIB) *CYP 

F6=2 . /A/B/ (SK2-4 . *PIA) *SD 

F7=2 . /A/B/ (SK2-9 . *PIB) *COS (3 . *PI*YP/B) 

C WRITE(6,7)F1,F2,F3,F4 

7 FORMAT ( 'Fl=' ,F10.5,5X, 'F2=' ,F10.5,5X, 'F3=' ,F10.5,5X, ' 
/F4=' ,F10.5) 

C WRITE(6, 8)F5,F6,F7 

8 FORMAT (' F5=' ,F10.5,5X, 'F6=' ,F10.5,5X, 'F7=' ,F10.5) 

DO 200 1=1, 181 

THDD=FLOAT(I-l) /2 . -45. 

THD=ABS (THDD) 

THR=THD/DPR 

PHR=PHD/DPR 

IF (THDD . LT . 0 . ) PHR=PHR+PI 
XS=RF*SIN(THR) *COS (PHR) 

YS=RF*SIN(THR) *SIN(PHR) 

ZS=RF*COS (THR) 

FX=CJJ 

FY=CJJ 

FZ=CJJ 

DO 104 14=1, N 
EX=CJJ 
EY=CJJ 
EZ=CJJ 

DO 102 12=1,4 
XA= (B+T) /2 . 

IF(I2 .EQ.2)XA=-XA 
IF (12 . EQ . 3 ) XA= (A+T) /2 . 

IF(I2.EQ.4)XA=- (A+T)/2. 

X=XS-XE(I4) -XA 
Y=YS-YE(I4) 

Z=ZS 

IF(I2 .LT.3)X=XS-XE(I4) 

IF(I2 .LT.3)Y=YS-YE(I4) -XA 
RR=SQRT (X*X+Y*Y+Z*Z) 

XU=X/RR 

YU=Y/RR 

ZU=Z/RR 

IF(ABS(XU) .LT.l.E-8)XU=l.E-8 
IF(ABS(YU) .LT.l.E-8)YU=l.E-8 
PHPR=ATAN2 (YU, XU) 

IF (PHPR . LT . 0 . ) PHPR=PHPR+TPI 
IF ( 12 . LT . 3 ) PHPR=PHPR- PI /2 . 

IF ( PHPR . LT . 0 . ) PHPR=PHPR+TPI 
THPR=ACOS(ZU) 
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CALL SL0T4 (XO , YO , ZO , XI , Y1 , Z1 , X2 , Y2 , Z2 , X3 , Y3 , Z3 ) 

IF(IANG(I4) .EQ.270)GO TO 50 

IF(IANG(I4) .EQ.90)GO TO 60 

IF(IANG(I4) .EQ.180)GO TO 70 

IF(I2.EQ.2)GO TO 41 

IF (12. EQ. 3) GO TO 42 

IF(I2.EQ.4)GO TO 43 

EX=EX- (F1-F2+F5-F7) *Y0- (F4-F3+F6) *Y2 
EY=EY+ (F1-F2+F5-F7) *X0+ (F4-F3+F6) *X2 
EZ=EZ+ (F1-F2+F5-F7) *Z0+ (F4-F3+F6) *Z2 
GO TO 102 

41 EX=EX+ (F2+F1+F5+F7) *Y0+ (F3+F4+F6 ) *Y2 
EY=EY- (F2+F1+F5+F7) *X0- (F3+F4+F6) *X2 
EZ=EZ- (F2+F1+F5+F7) *Z0- (F3+F4+F6) *Z2 
GO TO 102 

42 EX=EX+ (F1-F6) *X0+ (F2-F3) *X1+ (F4-F5) *X2-F7*X3 
EY=EY+ (F1-F6) *Y0+ (F2-F3) *Y1+ (F4-F5) *Y2-F7*Y3 
EZ=EZ+ (F1-F6) *Z0+ (F2-F3) *Z1+ (F4-F5) *Z2-F7*Z3 
GO TO 102 

43 EX=EX- (F1-F6) *X0+ (F3-F2) *X1- (F4-F5) *X2+F7*X3 
EY=EY- (F1-F6) *Y0+ (F3-F2) *Y1- (F4-F5 ) *Y2+F7*Y3 
EZ=EZ- (F1-F6) *Z0+ (F3-F2) *Z1- (F4-F5) *Z2+F7*Z3 
GO TO 102 

50 IF(I2 .EQ.2)GO TO 51 
IF(I2.EQ.3)GO TO 52 
IF(I2 .EQ.4)GO TO 53 

EX=EX- ( (F1-F6) *Y0+ (F3-F2) *Y1+ (F4-F5) *Y2+F7*Y3) 
EY=EY+( (F1-F6) *X0+ (F3-F2) *X1+ (F4-F5) *X2+F7*X3) 
EZ=EZ+( (F1-F6) *Z0+(F3-F2) *Z1+(F4-F5) *Z2+F7*Z3) 
GO TO 102 

51 EX=EX- (- (F1-F6) *Y0+ (F2-F3) *Y1- (F4-F5) *Y2-F7*Y3) 
EY=EY+ {- (F1-F6) *X0+ (F2-F3) *X1- (F4-F5) *X2-F7*X3) 
EZ=EZ+ (- (F1-F6) *Z0+ (F2-F3) *Z1- (F4-F5) *Z2-F7*Z3) 
GO TO 102 

52 EX=EX+ ( (F1-F2+F5-F7) *X0+ (F4-F3+F6) *X2) 

EY=EY+ ( (F1-F2+F5-F7) *Y0+ (F4-F3+F6) *Y2) 

EZ=EZ+ ( (F1-F2+F5-F7) *Z0+ (F4-F3+F6) *Z2) 

GO TO 102 

53 EX=EX+ ( (-F2-F1-F5-F7) *X0+ (-F3-F4-F6) *X2) 

EY=EY+ ( (-F2-F1-F5-F7) *Y0+ (-F3-F4-F6) *Y2) 

EZ=EZ+ ( (-F2-F1-F5-F7) *Z0+ (-F3-F4-F6) *Z2) 

GO TO 102 

60 IF(I2 .EQ.2)GO TO 61 
IF(I2 .EQ.3)GO TO 62 
IF(I2 .EQ.4)GO TO 63 

EX=EX- { (F1-F6) *Y0+ (F2-F3) *Y1+ (F4-F5) *Y2-F7*Y3) 
EY=EY+ (F1-F6) *X0+ (F2-F3) *X1+ (F4-F5) *X2-F7*X3 
EZ=EZ+ (F1-F6) *Z0+ (F2-F3) *Z1+ (F4-F5) *Z2-F7*Z3 
GO TO 102 

61 EX=EX+ (F1-F6) *Y0+ (F2-F3) *Y1+ (F4-F5) *Y2-F7*Y3 
EY=EY- ( (F1-F6) *X0+ (F2-F3 ) *X1+ (F4-F5) *X2-F7*X3) 
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EZ=EZ- ( (F1-F6) *Z0+ (F2-F3) *Z1+ (F4-F5) *Z2-F7*Z3) 

GO TO 102 

62 EX=EX+ (F1+F2+F5+F7) *X0+ (F3+F4+F6) *X2 
EY=EY+ (F1+F2+F5+F7) *Y0+ (F3+F4+F6) *Y2 
EZ=EZ+ (F1+F2+F5+F7) *Z0+ (F3+F4+F6) *Z2 
GO TO 102 

63 EX=EX+ (F2-F1-F5+F7) *X0+ (F3-F4-F6) *X2 
EY=EY+ (F2-F1-F5+F7) *Y0+ (F3-F4-F6) *Y2 
EZ=EZ+ (F2-F1-F5+F7) *Z0+ (F3-F4-F6) *Z2 
GO TO 102 

70 IF(I2 .EQ.2)GO TO 71 
IF(I2 .EQ.3)GO TO 72 
IF (12. EQ. 4) GO TO 73 

EX=EX- ( (F1+F2+F5+F7) *Y0+ (F3+F4+F6) *Y2) 

EY=EY+ (F1+F2+F5+F7) *X0+ (F3+F4+F6 ) *X2 
EZ=EZ+ (F1+F2+F5+F7) *Z0+ (F3+F4+F6) *Z2 
GO TO 102 

71 EX=EX- ( (F2-F1-F5+F7) *Y0+ (F3-F4-F6) *Y2) 

EY=EY+ (F2-F1-F5+F7) *X0+ (F3-F4-F6) *X2 
EZ=EZ+(F2-F1-F5+F7) *Z0+ (F3-F4-F6 ) *Z2 
GO TO 102 

72 EX=EX+ (F1-F6) *X0+ (F3-F2) *X1+ (F4-F5) *X2+F7*X3 
EY=EY+ (F1-F6) *Y0+ (F3-F2) *Y1+ (F4-F5) *Y2+F7*Y3 
EZ=EZ+ (F1-F6) *Z0+ (F3-F2) *Z1+ (F4-F5) *Z2+F7*Z3 
GO TO 102 

73 EX=EX+ (F6-F1) *X0+ (F2-F3) *X1+ (F5-F4) *X2-F7*X3 
EY=EY+ (F6-F1) *Y0+ (F2-F3) *Y1+ (F5-F4) *Y2-F7*Y3 
EZ=EZ+ (F6-F1) *Z0+ (F2-F3) *Z1+ (F5-F4) *Z2-F7*Z3 

102 CONTINUE 

FX=FX+EX*AMP(I4) *CEXP(-CJ*PHAS(I4) ) 

FY=FY+EY*AMP(I4) *CEXP(-CJ*PHAS(I4) ) 
FZ=FZ+EZ*AMP(I4)*CEXP(-CJ*PHAS(I4) ) 

104 CONTINUE 

ETH=FX*COS(PHR) *COS (THR) +FY*SIN (PHR) *COS (THR) -FZ*SIN 
/ (THR) 

EPH=-FX*SIN(PHR)+FY*COS (PHR) 

C ER=(ETH+CJ*EPH)/SQRT(2.) 

C EL= (ETH-CJ*EPH) /SQRT (2 . ) 

ER=ETH 

EL=EPH 

AER=CABS(ER) 

AEL=CABS(EL) 

IF(AER.LT.1.E-15)AER=1.E-15 
IF(AEL.LT.1.E-15)AEL=1.E-15 
IF(AER.LT.1.E-8)ER=CMPLX(1. ,1.) 

IF (AEL. LT . 1 . E-8) EL=CMPLX (1 . , 1 . ) 

DBR ( I ) =2 0 . *ALOGl 0 ( AER) 

DEL ( I ) =20 . *ALOG10 (AEL) 

PHASR ( I ) =ATAN2 (AIMAG (ER) , REAL (ER) ) *DPR 
PHASL (I ) =ATAN2 (AIMAG (EL) , REAL (EL) ) *DPR 
XAX(I)=THDD 



I 
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C WRITE(6, *)XAX,DBR,DBL 

200 CONTINUE 

YMAXL=-1000. 

YMAXR=-1000. 

DO 22 1=1,181 

IF(YMAXR.LT.DBR(I) ) YMAXR=DBR ( I ) 
IF(YMAXL.LT.DBL(I) ) YMAXL=DBL ( I ) 

22 CONTINUE 

YMAX=YMAXR 

IF ( YMAXR . LT . YMAXL) YMAX=YMAXL 
DO 21 1=1,181 
DBR(I)=DBR(I) -YMAX-2. 

DBL(I) =DBL(I) -YMAX-2 . 

IF(DBR(I) .LT. -40. )DBR(I)=-40. 

IF(DBLd) .LT.-40.)DBL(I)=-40. 

21 CONTINUE 

CALL BGNPLT 

CALL PLFORM( 'LINLIN' ,7.43,4.6) 

XX(1)=-180. 

XX(2)=180. 

YY(l)=-40. 

YY(2)=0. 

CALL PLSCAL(XX, 2, 080908, YY, 2, 040504) 

CALL PLGRAF( 'LP MICROSTRIP ARRAY' , 'ANGLE' , 'DB' ) 
CALL PLCURV(XAX,DBR,181,0,0) 

CALL PLNUP 
CALL PLNTYP(5) 

CALL PLNDN(0. ,0.) 

CALL PLCURV(XAX,DBL, 181, 0, 0) 

CALL ENDPLT 

J=181 
T=-181 

DO 211 1=182,362 
XAX(I)=(I-l)/2.-45. 

IF(XAXd) .GT.90. )XAXd)=XAX(I) -180. 

DBR(I) =DBR(J) 

J=J-1 
211 CONTINUE 

OPEN (38, F ILE=' RESULTS. PRN ' ) 

WRITE (6, 699) 

699 FORMAT ( ■ WRITING RESULTS TO FILE. . . ' ) 

700 F0RMAT(F8.3,1X,F8.3,1X) 

DO 701 1=1,361 
WRITE(38, *)XAX(I) ,DBR(I) 
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CONTINUE 



WRITE (6, *) ' DONE. . . ' 

1000 CLOSE (33) 

CLOSE (35) 

CLOSE (38) 

STOP 

END 

SUBROUTINE SLOT4 (EX, EY, EZ , EXl , EYl , EZl , EX2 , 
/EY2 , EZ2 , EX3 , EY3 , EZ3 ) 

C. . .RADIATION FROM DIFFERENT MODE SLOTS 

COMPLEX C J , FAC , EPH , ETH , EX , EY , EZ , C JJ 

COMPLEX EXl , EYl , EZl , EX2 , EY2 , EZ2 , EX3 , EY3 , EZ3 

COMMON/DDl/EPS,T,AA,BB, 12, IP,XS,YS, ZS,CJJ 

COMMON/DD2/CJ, PI , TPI , DPR 

COMMON/DD4/THR, PHR, R 

FAC=CEXP ( -CJ*TPI*R) /R 

A=T/2 . 

B= (AA+T) /2 . 

IF ( 12 . GT . 2 ) B= (BB+T) /2 . 

CPH=COS (PHR) 

SPH=SIN(PHR) 

CTH=COS (THR) 

STH=SIN(THR) 

ARG1=TPI*A*CPH*STH*SQRT (EPS) 

ARG2 =TP I * B * S PH* STH 
IF(ABS(ARG1) .LT.l.E-4)GO TO 11 
F1=SIN (ARGl) /ARGl 
GO TO 12 

11 Fl=l. 

12 IF(ABS(ARG2) .LT.l.E-4)GO TO 13 
F2=SIN (ARG2 ) /ARG2 

GO TO 14 

13 F2=l. 

14 IF(ABS(ABS(ARG2)-Pl/2.) .LT.l.E-4)GO TO 15 
F3=ARG2*COS (ARG2) / (ARG2*ARG2-PI*Pl/4 . ) 

GO TO 16 

15 F3=-1./PI*ARG2 

16 F4=COS(ARG2)/(ARG2*ARG2-PI*PI) 
IF(ABS(ABS(ARG2) -1.5*PI) .LT.l.E-4)GO TO 17 
F5=ARG2*COS (ARG2) / (ARG2*ARG2-PI*PI*9 . /4 . ) 

GO TO 18 

17 F5=-l./2. 

18 EPH=-CJ*F1*F2*SPH*CTH*FAC 
ETH=CJ*F1*F2*CPH*FAC 
EX=ETH*CTH*CPH-EPH*SPH 
EY=ETH*CTH*SPH+EPH*CPH 
EZ=-ETH*STH 
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EPH=-F1*F3*SPH*CTH*FAC 

ETH=F1*F3*CPH*FAC 

EX1=ETH*CTH*CPH-EPH*SPH 

EY1=ETH*CTH*SPH+EPH*CPH 

EZ1=-ETH*STH 

EPH=CJ*PI*F1*F4*SPH*CTH*FAC 
ETH=-CJ*PI*F1*F4*CPH*FAC 
EX2=ETH*CTH*CPH-EPH*SPH 
EY2=ETH*CTH*SPH+EPH*CPH 
EZ2=-ETH*STH 
EPH=F1*F5*SPH*CTH*FAC 
ETH= - FI * F5 *CPH* FAC 
EX3=ETH*CTH*CPH-EPH*SPH 
EY3=ETH*CTH*SPH+EPH*CPH 
EZ3=-ETH*STH 
RETURN 
END 



III. GEOM.FOR 

C. . .ROUTINE TO CREATE GEOMETRY FOR MI CARY. FOR 

DIMENSION X(4) ,Y(9) ,SMAGD{9,4) , SPHASED{9, 4) 
REAL L,W, AMP, PHASE, LAMBDA, FREQ, SPACE 

OPEN ( 5 , FILE= ' GEOMETRY ' ) 

OPEN(33,FILE='MAG' ) 

OPEN (35 , FILE= ' PHASE ' ) 

C. . .READ IN MAGNITUDE AND PHASE DATA FROM FILE 



NUMX=4 

NUMY=9 

NUMTOT=NUMX*NUMY 
DO 4 I=1,NUMY 

READ(33, *)SMAGD(I, 1) ,SMAGD(I,2) ,SMAGD(I,3) ,SMAGD(I,4) 
READ(35, *)SPHASED(I, 1) , SPHASED(I, 2) , SPHASED(I, 3) , 
/SPHASED (1,4) 

WRITE(6,*)SMAGD(I,1) ,SMAGD(I,2) ,SMAGD(I,3) ,SMAGD(I,4) 
WRITE(6,*) SPHASED (1,1) , SPHASED ( I , 2 ) , SPHASED ( I , 3 ) , 
/SPHASED (1,4) 

4 CONTINUE 

12 FORMAT(F10.5,F10.5,F10.5,F10.5) 

C . . . FREQUENCY IN HZ 
FREQ=1618.25E6 
LAMBDA=3 . E8/FREQ 
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C. . .DIMENSIONS IN CM 
L=4 .44 
W=5.9 
AMP=1 . 0 
PHASE= . 0 
ORIENT=0. 

SPACE= . 55*LAMBDA 

C . . . CONVERT TO WAVELENGTHS 
L=L/100. 

L=L/LAMBDA 

W=W/100. 

W=W/LAMBDA 

DO 10 I=1,NUMX 
X(I)=(I-1)*SPACE+(I-1) *W 
10 CONTINUE 

DO 15 I=1,NUMY 
Y(I)=(I-1) *SPACE+(I-1) *L 
15 CONTINUE 

NUMEL=NUMX*NUMY 
WRITE(5, 98)NUMEL 

98 FORMAT (13) 

K=0 

N=0 

ORIENT=0 
DO 20 1=1,9 
DO 22 J=l,4 

K=K+1 

WRITE(5,99)X(J) ,Y(I) ,SMAGD(I,J) , SPHASED (I , J) ,N 

22 CONTINUE 
20 CONTINUE 

99 FORMAT (F8 . 3 , IX, F8 . 3 , IX, F8 . 3 , IX, F8 . 3 , IX, II) 

STOP 

END 
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APPENDIX G 



I . RADIATION PATTERNS 



i 




Figure G-1. Beam 1 Radiation pattern (<|)=90®) 




Figure G-2. Beam 2 Radiation Pattern ((j)=90°) 
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Figure G-3. Beam 3 Radiation Pattern (<|)=90®) 







Figure G-4. Beam 4 Radiation Pattern (<|)=90®) 
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Figure G-5. Beam 5 Radiation Pattern ((|)=90®) 
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II 



PATTERN SEPARATION 



0 




Figure G-7. Simultaneous Radiation 
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Figure G-8. Simultaneous Radiation 
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Figure 6-9. Simultaneous Radiation 



III. ANTENNA FOOTPRINT 



Relative dB Contours 
(Beam scanned 45 degrees from nadir) 




Figure G-10. Beam 1 Relative dB Gain Footprint 
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Relative dB Contours 
(Beam scanned 27 degrees from nadir) 




Left Right 

Cross Track 



Figure G-11. Beam 2 Relative c3B Gain Footprint 



Relative dB Contours 
(Beam scanned 10 degrees left of nadir) 




Left Right 

Cross Track 



Figure G-12. Beam 3 Relative dB Gain Footprint 
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Relative dB Contours 
(Beam scanned 10 degrees right of nadir) 




Figure G-13. Beam 4 Relative dB Gain Footprint 



Relative dB Contours 
(Beam scanned 27 degrees from nadir) 

9 0_, 




Left Right 

Cross Track 



Figure G-14. Beam 5 Relative dB Gain Footprint 
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Relative dB Contours 
(Beam scanned 45 degrees from nadir) 




Left Right 

Cross Track 



Figure G-15. 



Beam 6 Relative dB Gain Fcxstprint 
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